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Fig. 2. — (a) The fraction of the critical density in baryons (hydrogen and helium) associated 
with damped Lya systems at z = 3.2. The thick solid curves (top for Q v =0. 2, bottom for Q v =0. 3) 
assume the hydrogen column density threshold of Lanzetta et al. (1993). The top two dotted curves 
show Q g predicted by the Q v =0. 2 model, assuming all hydrogen in Denmax halos with 6p/p > 200 
(bottom) and > 50 (top) contributes to Lya absorption; the bottom two are for Q v =0. 3. The 
shaded band shows the observed Q g from Lanzetta et al. (b) The column density distribution 
predicted by the Q v =0. 2 model (thick line) vs. data points from Lanzetta et al. for 2.5 < z < 3.5. 
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Fig. 1. — The cumulative comoving number density of dark halos as a function of z for the two 
CDM+HDM models. In each panel, the five solid curves show the Press-Schechter fits with a tophat 
filter and 6 C (indicated) above five mass thresholds: M = 10 13 , 10 12 ' 5 , 10 12 , 10 11 ' 5 , and lO n M 
(from bottom up). The filled circles joined by dotted lines represent the simulated halos above the 
same mass thresholds. The Press-Schechter function for the CDM model (6 C = 1.8 for both plots) 
is also shown for M > 10 11 (top dashed) and > 1O 13 M0 (bottom dashed). The five filled squares 
with error bars represent the comoving space density of quasars (for qo = 0.5 and Q = 1) brighter 
than M B = -26 from Schmidt et al. (1991, 1994). 
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either case our approach is more reliable than estimates based on analytic approximations such as 
the Press-Schechter formula combined with assumptions such as universal velocity-dispersion/halo 
mass relations. Our procedure likely overestimates Q g because we assume that all hydrogen 
contributes to the absorption (i.e., there is no ionization and no dark baryons). Star formation 
and ionization by shocks and ultraviolet background radiation will further reduce our estimate of 
Q g in the models. 

Figure 2a shows the striking discrepancy between the Q v = 0.3 model and the observation at 
z ~ 3.2: only 0.007% of the total mass density in this model is in the form of high column- density 
baryonic clouds compared to 0.058% from observations. For Q v = 0.2, the fraction rises to 
~ 0.02%, still a factor ~ 3 below the value needed to match the observation, even with our 
conservative assumptions. The column density distribution predicted by the model is also too 
steep (Fig. 2b). Ionization and gaseous dissipation will decrease the slope of f(N) but ionization 
will also decrease Q g , making the discrepancy worse. If the power spectrum normalization is 
increased to (Jrms-PS = 21 /xK (2.3<r above the estimate of Bennett et al. 1994), Q g is increased by 
a factor ~ 1.5, still not enough for Q v = 0.3 to match the observations. We thus conclude that if 
damped Lya systems are indeed embedded in dark halos with a primordial baryon fraction, then 
their high abundance constrains the neutrino mass density in CDM+HDM models to be Q v £, 0.2, 
providing an upper bound of 4.7 eV on the most massive light neutrino. 

Klypin et al. (1994) have reached a different conclusion. However, their Figure 2 indicates 
that even all the galaxies with total mass exceeding 1O 1O M0 do not provide enough Lya absorption 
at z ~ 3.2 for Q v > 0.25. Our numerical calculations (Fig. 2) suggest that the damped Lya 

f . 



absorbers in CDM+HDM models with Q v > 0.2 come from halos more massive than 5 X 10 n M(: 



Our results place very strong constraints on the CDM+HDM models. While the problems of 
high-redshift galaxy formation can be ameliorated by decreasing the HDM fraction further, this 
has the effect of increasing the masses and relative velocities of galaxies at z = 0. Klypin et al. 
(1993) and Nolthenius, Klypin, & Primack (1994) concluded that the Q v = 0.3 model does not 
produce excessive small-scale power at z = 0, but questions remain for models with smaller Q v . 
In a later paper we shall present our own results on the properties of Q v = 0.2 and 0.3 models at 
z = 0. 

We are grateful to Maarten Schmidt for providing the current grism survey data prior to 
publication. We thank Sterl Phinney and Chuck Steidel for helpful comments on the manuscript. 
C.-P. Ma acknowledges support from a PMA Division Fellowship at Caltech. E. B. acknowledges 
support from NSF grants AST90-01762 and ASC93-18185. Supercomputer time was provided by 
the National Center for Supercomputing Applications. 



mass M > 7 X 10 8 M & . The black hole forms presumably by collapse of baryons, which in our 
models make up only 5% of the total mass. Thus, assuming 1% of the baryons go into the hole 
(a reasonable but uncertain estimate), the minimum galactic mass required to produce such a 
hole is 1.4 X 10 12 M & . In order to account for the observed numbers of quasars, the Q, v = 0.3 
model requires host galaxy masses at least 3 times smaller from above. We have assumed that 
the time lag between galactic halo formation and black hole formation is negligible; relaxing this 
assumption only increases the difficulty (Turner 1991). A constant quasar space density from 
z ~ 2 to 4 (Irwin et al. 1991) would also enlarge the discrepancy. By contrast, the Q v = 0.2 
model appears compatible with the high-redshift quasar counts. However, we cannot completely 
eliminate the Q v = 0.3 model based on quasar abundances because it is possible, in principle, 
that black holes form quickly and accrete more than 1% of the baryonic mass in galaxies hosting 
quasars. 

An independent test comes from the amount of neutral gas in galaxies associated with 
damped Lya systems. Lanzetta (1993) and Lanzetta, Wolfe, & Turnshek (1993) estimated that 
the fraction of the critical density in the dense neutral gas associated with damped Lya systems 
is (for qo = 0.5) £l g = (5.8 ± 1.2) X 10~ 3 h$ at 3.0 < z < 3.5, a value that can be compared with 
predictions of models. 

To calculate the amount of dark mass in our simulations that is associated with the observed 
damped Lya systems, we first compute for each of our Denmax halos (without any overdensity 
cut) the mass column density as a function of projected radius from the halo center. We obtain 
the corresponding hydrogen column density JV(HI) by multiplying the dark mass column density 
by f2{,(/xmH) _1 (where mH=hydrogen mass; /x=1.33 to take into account the 25% of helium by 
mass) under the assumption that the gas and dark matter are distributed similarly in the outer 
parts of the halos. To compare fairly with observations, we then compute for each halo only the 
baryonic mass (including both hydrogen and helium) in the cylinder whose JV(HI) exceeds the 
corresponding completeness threshold 2 X 10 20 cm -2 in the survey of Lanzetta et al. (1993). Figure 
2a shows Q g computed in this way at z = 3.2 for halos more massive than M. For comparison 
we show Q g computed assuming that all the baryons in halos overdense by factors of 50 and 
200 contribute to damped Lya absorption. Figure 2b compares the column density distribution 
function f(N) [related to Q g by Q g = i?oc _1 /xmH/0~ ;t f^ . N f(N)dN] observed by Lanzetta et al. 
with the predictions of the Q v = 0.2 model. 



The only assumption behind our calculation is that the damped Lya systems are embedded 
in dark halos where the gas and dark matter are distributed alike with a constant ratio given by 
primordial nucleosynthesis. Because the corresponding halo radii exceed 40 kpc for M > 1O 12 M0, 
dissipative effects should be small so that our assumption is reasonable. Alternatively, if all 
baryons in dark halos dissipate to high column density, the reader may refer to the dotted lines 
in Figure 2a. Note that the thick lines in Figure 2a in fact give higher Q g for M <; IO^Mq than 
the dotted lines for virialized halos. This is because we conservatively used all Denmax halos (and 
not just those with 8pj p > 200) for the thick lines when we computed the column densities. In 



Schechter 1974) for the number density of virialized mass lumps: 
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where po(z) denotes the background mass density of the universe, 8 C is a free parameter 
characterizing the linear overdensity at the onset of gravitational collapse, and <r 2 (M, z) is the 
variance of the density field at redshift z, computed using the linear power spectrum P(k,z) 
and a smoothing window function W(kR). The power spectrum used here is taken to be the 
density- averaged P = [£l v P v + (f2 c + f2{,)P c ] 2 , where P v and P c are computed numerically from 
the coupled, linearized Boltzmann, fluid, and Einstein equations that govern the evolution of the 
CDM and HDM density fields (e.g., Ma & Bertschinger 1994b). 

We find that the best-fit parameter for the Q v = 0.3 model is 8 C ~ 1.8 for a spherical tophat 
window function instead of the canonical value 1.68 (Fig. 1). This difference in 8 C significantly 
changes the predicted abundance for rare objects at high redshift. (A Gaussian window function 
requires a smaller 6 C , in this case 8 C ~ 1.7. The tophat gives a somewhat better match in the 
shape.) The Q v = 0.2 model requires different 8 C for early and late epochs. At z ^ 3, 8 C ~ 1.8 
provides a reasonable fit, but this choice underpredicts the halo abundance by a factor of 3 at 
z = 5.3 and 6 at z = 7.4. At z > 3, 8 C = 1.68 provides a better fit (Fig. 1). This trend of favoring 
a smaller 8 C at earlier times was also observed in the pure CDM model by Gelb & Bertschinger 
(1994). Overall, however, the slope dlogN/dlogM is steeper in the Press-Schechter approximation 
than in our simulations. 



Observational Constraints 



The comoving number density of quasars brighter than Mb = — 26 from the Palomar grism 
survey (Schmidt et al. 1991) is plotted in Figure 1. The data point at z=2.2 is based on 
C IV-selected quasars; the points at higher z are from the recent analysis of Lya-selected quasars 
(Schmidt et al. 1994). The highest-redshift (z = 4.5) data point provides the most stringent 
constraint. For Q v = 0.3, no massive halos have formed at this redshift in our simulation box, 
so the theoretical predictions can be estimated only by extrapolating the best-fit Press-Schechter 
curves to our simulation results. If quasars are active for only a fraction / of the time, then the 
required number of host halos increases by a factor / _1 . Even in the most optimistic case (/ = 1), 
when 100% of the dark matter halos at high redshift host quasars at any one time, our numerical 
results show that every halo of mass greater than ~ 5 X 10 11 M & would have to host a luminous 
quasar for this model to match the observations. For Q v = 0.2, this mass rises to ~ 7 X 10 12 Mq. 

We can assess whether these requirements are plausible by assuming that a quasar corresponds 
to a massive black hole radiating at the Eddington limit. With a bolometric correction factor 
of 6.0 (Sanders et al. 1989; Haehnelt & Rees 1993), Mb < —26 corresponds to a black hole 



begun at a redshift z = 13.6 (Q v = 0.3) and z = 21.5 (Q v = 0.2), when the rms density fluctuation 
in the CDM component was about 0.2. 

The simulations of the two models were performed with identical parameters (aside from the 
neutrino masses) and the same realization of the Gaussian random field. The primeval fluctuations 
are assumed scale-invariant and isentropic, normalized to a COBE rms quadrupole moment 
<5rms-PS = 17.6 /xK (Bennett et al. 1994). The simulation volume is a cube of comoving length 
100 Mpc and the interparticle force law is a Plummer law with a comoving softening length of 50 
kpc (or a proper softening length of 10 h$ kpc at z = 4). We used 128 3 simulation particles to 
represent the CDM and 10 times as many HDM particles for an accurate sampling of the neutrino 
phase space. For the Q, v = 0.2 case the particle masses are 2.5 X 10 10 M & for CDM and 6.6 X 10 8 
M & for HDM. The simulations evolved to z = 1.5 required a total of 1300 CPU hours on a Convex 
C-3880. 

The dark matter halos in our simulations are identified with an improved version of the 
Denmax algorithm described in Gelb & Bertschinger (1994). At the redshift of interest, Denmax 
finds all the particles within closed density contours around each density peak by moving the 
particles along the gradient of the mass density field (defined using a Gaussian of comoving radius 
100 kpc) until they settle into the peaks. This method identifies distinct density concentrations 
more reliably than the standard friends-of-friends algorithm. After identifying the particles 
associated with each peak we compute the energy of the particles in each halo and remove the 
unbound ones. CDM and HDM particles are treated identically. 

Because the task of Denmax is to associate particles with density peaks regardless of the 
value of the density maximum, some of the halos identified by Denmax at high redshifts have not 
yet collapsed and virialized. We therefore impose an overdensity constraint of 6p/ p > 200 to select 
virialized halos. The effect of this cut is to reduce the mass of all halos and to eliminate small 
density enhancements that have not collapsed. The halo abundance is sensitive to the overdensity 
constraint. We find up to 5 times more halos when we apply the less stringent constraint 6p/ p > 50 
used by Klypin et al. (1993; see Fig. 2a below). The higher overdensity constraint is appropriate 
for virialized halos. 

The comoving number density of Denmax halos more massive than M, N(M,z), is shown in 
Figure 1 for the two models. It is a very sensitive function of Q v and Q c at high redshifts. The 
results at z > 4 are not plotted for Q v = 0.3 because the slower growth in this model produces 
little structure at this early time, and it is difficult to identify halos reliably. In contrast, virialized 
halos with M ~ lO n ' 5 M0 are present as early as z ~ 7.4 in our 100 3 Mpc 3 simulation box for 
Q v = 0.2. In fact, the halo abundance at z = 7.4 in the Q v = 0.2 model is nearly identical to that 
at z = 4.0 in the Q v = 0.3 model. 

We can calibrate and check the validity of the Press-Schechter approximation (Press & 



of quasars at z > 2. The Palomar CCD grism survey (Schmidt, Schneider, & Gunn 1991, 1994) 
has yielded a decline in the space density of quasars with Mb < —26 by a factor ~ 3 per unit 
redshift for z <; 2.5. Multicolor photographic surveys, on the other hand, have shown a slightly 
steeper decline in one case (Warren, Hewett, & Osmer 1994) and are consistent with a constant 
comoving density for Mb Ss —28 at 2 < z < 4 in another (Irwin, McMahon, & Hazard 1991). 

Damped Lya absorption systems offer an alternative measure of early galaxy formation. 
Whether these systems are associated with gas-rich dwarf galaxies or normal disk spirals is not 
absolutely clear. However, recent observations of QSO H0836+113 (Wolfe et al. 1992) have failed 
to confirm the detection of a narrow emission feature (Hunstead et al. 1990) which had been 
regarded as evidence for a star-forming dwarf. Wolfe et al. (1992) instead detected extended line 
emission from a large area of 96 X 44 h$ kpc 2 in the same system. Current surveys appear to 
favor the association of high column density Lya absorbers with normal disk galaxies instead of 
dwarfs (e.g., Steidel & Hamilton 1992, 1993). 

Assuming that the damped Lya systems are mainly dwarf galaxies, Klypin et al. (1994) 
concluded that the CDM+HDM models are viable, in contrast with the analytic results of Mo & 
Miralda-Escude (1994), Kauffmann & Chariot (1994), and Subramanian & Padmanabhan (1994). 
Instead of relying on analytic approximations, we will estimate directly from large cosmological 
JV-body simulations the fraction of the mass with hydrogen column density above the survey 
completeness threshold and compare with observations. We will also compare the abundance of 
dense dark matter halos capable of hosting quasars with the observed quasar abundance. 

2. Simulations and Halo Properties 



We have performed high-resolution CDM+HDM gravitational JV-body simulations using 
the adaptive particle-particle/particle-mesh code (P 3 M) 2 (Bertschinger & Gelb 1991; Gelb & 
Bertschinger 1994). In all cases Q v + Q c + Q^ = 1 and Ho = 50 km s _1 Mpc -1 . Two different 
neutrino masses were considered: m v = 7.0 eV (Q v = 0.3) and m v = 4.7 eV (Q v = 0.2). The 
baryon fraction used in computing the initial conditions was set to the value implied by primordial 
nucleosynthesis, Q^ = 0.05 (Walker et al. 1991). In the JV-body simulations the baryons were 
treated identically as CDM. 

The initial conditions for the hot and cold particles were generated using the Monte Carlo 
procedure described in Ma & Bertschinger (1994a). Our method differs from that used by other 
workers (Davis et al. 1992; Klypin et al. 1993; Cen & Ostriker 1993; Jing et al. 1994) in that 
the initial positions and momenta of our HDM particles were sampled from the full perturbed 
neutrino phase space distribution, while their HDM momenta were drawn from the Fermi-Dirac 
distribution which is the zeroth-order approximation to the full distribution function. We also 
took into account the wavenumber dependence of the suppression of the CDM velocities owing to 
neutrino damping in these models (Ma & Bertschinger 1994a). The nonlinear calculations were 



1. Introduction 

The cold+hot dark matter (CDM+HDM) model of structure formation has emerged as a 
promising alternative to the standard CDM model. Replacing some of the CDM with HDM 
(massive neutrinos) suppresses small-scale power for a fixed large-scale normalization such as that 
provided by the cosmic microwave background anisotropy. This suppression helps to alleviate the 
problem of excessive small-scale power in the CDM model (e.g., Gelb & Bertschinger 1994) but it 
also delays the epoch of galaxy formation. Thus, the abundance of high-redshift quasars and the 
amount of mass in dense gas clouds can be used to constrain the amount of HDM present. 

Early works based on linear calculations (Schaefer, Shah, & Stecker 1989; van Dalen 
& Schaefer 1992; Taylor & Rowan- Robinson 1992; Holtzman & Primack 1993) and JV-body 
simulations (Davis, Summers, & Schlegel 1992; Klypin et al. 1993; Jing et al. 1994) have shown 
general agreement between observations at redshifts z ~ and the Q v = 0.3 flat CDM+HDM 
model. Some of these workers have expressed concern about late galaxy formation, but the lack 
of adequate dynamic range in mass and length in the fully nonlinear regime prevented them from 
reaching firm conclusions. 

Analytical approximations have been used by many workers to avoid the limited dynamic 
range of numerical simulations. The abundance of bound objects at high redshifts has been 
estimated using the Press-Schechter (Kashlinsky 1993; Haehnelt 1993) and the Gaussian peak 
(Nusser & Silk 1993) formalism in the hope of constraining various dark-matter models using 
quasar statistics. Similar calculations have been made comparing analytic models with the damped 
Lya absorption systems (Mo & Miralda-Escude 1994; Kauffmann & Chariot 1994; Subramanian & 
Padmanabhan 1994). However, the accuracy of these approximations in the appropriate regime is 
uncertain. The peak model for halo formation, for example, has been criticized recently by Katz, 
Quinn, & Gelb (1993). As we will show, the Press-Schechter formula also can lead to inaccurate 
results for rare objects if not properly calibrated. 

To be sure, uncertainties remain in various aspects of both the theories and observations of 
quasars and Lya systems. Our understanding of the relation between the observed quasar activity 
on parsec scales and the host galaxies of size ~ 10 kpc is incomplete. Quasar formation itself 
cannot be addressed directly by any current cosmological simulations, but progress has been made 
in modeling the formation and fueling of black holes in the host galaxies of quasars. Efstathiou & 
Rees (1988) observed that the epoch of galaxy formation in the biased (b ~ 2.5) CDM model is 
consistent with the observed luminous quasar density between z = 2 and 4, but a sharp cut-off 
at z > 5 is predicted. Turner (1991) argued that there is a considerable time lag between the 
formation of virialized halos and the birth of high-redshift quasars, pushing the halo formation 
epoch to z > 5 or possibly z > 10. A contrary view was recently expressed by Haehnelt & Rees 
(1993), who argued for a short delay between halo collapse and formation of a massive black hole. 

Uncertainties also exist on the observational front in estimating the comoving number density 
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The abundance of galactic systems at high redshifts can impose a strong constraint 
CN on the cold+hot dark matter (CDM+HDM) models. The hot component reduces 

the excessive small-scale power in the COBE-normalized CDM model but also 
delays the epoch of galaxy formation. We present results from the first numerical 



(^ simulations that have enough dynamic range to address accurately the issue of 



high-redshift halo abundances in CDM+HDM models. Equivalent high-resolution 



■^- particle-particle/particle-mesh JV-body simulations are performed for spatially flat 

models with 0„ = 0.3 and 0.2 (with H = 50 km s" 1 Mpc" 1 and Q b = 0.05). We 

'■3 study the constraints placed on the models by the high-redshift quasar space density 

i and by the mass fraction in neutral dense gas associated with damped Lya systems. 

5^ We find that even with optimistic assumptions, the much-studied Q v = 0.3 model does 

not produce enough massive halos to account for the observed abundance of quasars 
at z > 4. The model passes this test if Q v is decreased to 0.2. Both models do not 
produce enough high column- density halos to account for the amount of gas in damped 
Lya systems at z <; 3: the Q v = 0.3 model falls short by a factor ~80; the Q v = 0.2 
model by a factor ~3. We conclude that only CDM+HDM models with Q v <^ 0.2 can 
match observations at high redshift, implying an upper bound of 4.7 eV on the most 
massive light neutrino (presumably the r). 
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